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ABSTRACT 

The next generation of astrometric instruments will reach accuracies deserving new 
treatments. In order to get astrometric parameters achieving the precision permitted 
by the measurements, it will be necessary to take into account effects that were ne- 
glected until the present time. Two effects concerning the orbital elements of binary 
stars are considered hereafter: the former is the local perspective (LP) effect, which 
is due to the variation of the distance and of the orientation of the orbital plane dur- 
ing the observation time span; the equations describing this effect are derived for the 
first time. The latter effect is the light-travel time (LTT), which is also related to the 
orientation of the orbital plane, and which is as efficient as the preceding one. Taking 
these effects into account would allow to find the ascending nodes of the orbits, and 
lead to orbital elements more accurate than when they are ignored. 

It is derived from simulations that, at a distance of 5 pc, and assuming velocities 
typical of Pop. I stars, the position of the right ascending node could be derived for 
a few simulated unresolved binaries when the astrometric measurements have errors 
around 1 nas. For the resolved brown dwarf binary 2MASS J07464256 +2000321, it 
appears that ignoring the LP effect would result in underestimating the masses of 
the components by 14 per cent of the errors as soon as the astrometric errors are 
around 20 p.as for each measurement. However, a 'degenerate LP solution', taking 
into account the variation of the semi-major axis when the distance is varying, should 
provide reliable masses when the measurement errors are larger than 1 or 2 |j.as. A 
few binaries in the Gaia program could deserve a degenerate LP solution, whereas a 
the complete LP+LTT solution could be justified for resolved binaries observed with 
SIM. 

Key words: astrometry - methods: analytical - stars: binaries: general - stars: 
individual: 2MASS J07464256 +2000321 - stars: fundamental parameters 



1 INTRODUCTION 

Thanks to the progresses in spatial astrometry, the errors of 
measurements of the positions of stars are gradually decreas- 
ing: the errors of the along-scan measure ments were a few 
milli-arcseconds in the Hipparcos p roject {ESA 1997]), but 
the Gaia satellite JPerrvmanl 20051 ) is expected to achieve 
100 times better (|Mignardl 120031 1 . This limit will still be 
pushed beyond 1 (j.as when interferometry with large base- 
line will be feasible in space. It is expected that, with a 9-m 
baseline, SIM could alread y get an accuracy of only 0.6 [ias 
for a single measurement ( Unwin et al.|[200a ). 

These important improvements will not only lead to 
more accurate results. They make also necessary to take 
into account effects that were never considered in the past. 
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A first example was the correction for perspective accelera- 
tion, that was introduced for a few high-velocity single stars 
in the Hipparcos catalogue. The radial velo cities (RV) of 
these stars were used for that purpose, but iDravins et al.l 
l| 19991 ) ha ve shown that the RV may also be derived from 
this effect. IKlionerl i|2003l ) developed a relativistic model tak- 
ing into accou nt the light deflection by solar system objects. 
IKlionerl (|2000| ) and iMignardl (|2003h outlined that some ba- 
sic parameters as the trigonometric parallax, the proper 
motion and the RV deserve to be defined more precisely 
than they are at present; the d e finitio n of the RV was dis- 
cusse d by Lindegren fc Dravinl ([2003 ). Concerning binary 
stars, lAnglada-Escude fc Torral (|2006h have shown that the 
variations of the light-travel time (LTT) coming from the 
orbital motion may change the apparent position of the com- 
ponent. Since this effect is related to the radial displacement 
of the component, they inferred that it could be used to dis- 
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sipate the ambiguity of the position of the ascending node 
of the orbit. Another effect which could be used for node 
identification is the modification of the orientation of the 
orbit during the time span of the astrometric program. For 
sake of brevity, it is called 'local perspective' (LP) hereafter. 
In order to avoid confusion, the perspective effects concern- 
ing the barycentre of binary systems, but also single stars 
(i.e. perspective acceleration and change in the parallax), 
are called 'barycentric perspective' hereafter. 

In the course of a high-precision astrometric project, 
taking perspective and LTT effects into account could be 
necessary for two reasons. First of all, when the effects are 
relevant, the quality of the solution is improved. The as- 
trometric and orbital parameters are then more accurate, 
and, even more important, the solution itself will look more 
'acceptable'. The decision to consider a model as adequate 
is usually based on the goodness-of-fit (GOF) of the solu- 
tion, which may, for instance, follow the normal distribution 
N(Q, 1) when the model is the right one. Therefore, a large 
GOF leads to the rejection of the model, and to the search of 
another one. In the course of a large astrometric survey, like 
Hipparcos, objects which are not satisfactorily fitted even- 
tually got a 'stochastic' solution, which is of rather poor 
interest. Therefore, a binary model neglecting perspective 
or LTT effects could change actual binary stars in stochas- 
tic objects, losing the orbital elements. Aside from giving 
the right classification with accurate parameters, the sec- 
ond advantage of giving a precise solution is the possibility 
to get more informations about the object. For single stars, 
as mentioned above, taking into account perspective accel- 
eration can lead to an estimation of the RV. For binaries, a 
consequence of the additional effects could be the complete 
derivation of the orientation of the orbit in space. This is 
investigated hereafter. 

Calculations were performed in order to answer three 
questions: What is the smallest astrometric error still per- 
mitting to neglect LP and LTT without altering dramati- 
cally the GOF for unresolved binaries? What is the largest 
astrometric error permitting to derive the precise ascending 
node thanks to LP and LTT? What is the smallest astro- 
metric error still permitting to neglect LP and LTT without 
altering the masses of resolved binary components? The ba- 
sic equations of LP are derived in Section and it is shown 
how to select the right ascending node thanks to this effect. 
A method for selecting the ascending node from LTT is pre- 
sented in Section |31 In Section U it is verified that the rela- 
tivistic effects, which were neglected in the calculations, are 
much smaller than LP and LTT effects for nearby binaries 
with large semi-major axes. The parameters of simulated 
unresolved astrometric binaries were derived in Section [5] 
and the ascending nodes were searched. Resolved binaries 
are treated in Section HJ] Section [7] is the conclusion. 



2 THE LOCAL PERSPECTIVE EFFECT 
2.1 General transformations 

We want to derive the angular coordinates of binary com- 
ponents in a tangential reference frame moving with the 
barycentre. At the beginning, the components are located in 
the 'local equatorial reference system' of the object, which is 
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Figure 1. The first coordinate transformation in computing the 
local perspective effect: Rotation of angle of the local equatorial 
reference frame around the radial vector for aligning the rectan- 
gular coordinates with respect to the proper motion. The angle 8 
is the position angle of the proper motion; it is measured in the 
anticlockwise direction, but inside the celestial sphere. 
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Figure 2. Perspective transformation when the barycentre of the 
object is moving in the 3-D proper motion direction. The reference 
system is rotating with the a angle (er = p, t, or cr = p, t X (1 — 
p r t) when perspective acceleration is taken into account), and the 
angular dimensions are transformed by a factor d/d' = 1 — fi r t. 
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defined as a trihedron having the origin on the barycentre of 
the binary, and with the vectors k) oriented as follows: 
k is prolongating the line of sight, i is oriented toward east, 
and j is oriented to north, in order to complete the direct 
trihedron (figure [T]). 

The angular coordinates are lengths divided by the dis- 
tance from the Sun. At the epoch t = 0, the object is moving 
in space with the angular velocity ([i a *, fi$, fir). The third 
coordinate, fi r is the 'radial proper motion', which is always 
expressed in radians per time unit. It may be derived from 
the RV, v r , and from the parallax, zu, using the equation 
fi T = v r zu x (1 — v T /c), where v r is in AU per time unit and 
Tjo is in radians; c is the speed of light, in the same unit as 

Mr. 

We consider a point with a fixed position in an inertial 
reference frame attached to the barycentre (for instance: the 
periastron of an astrometric orbit). When the binary is mov- 
ing, the angular coordinates of this point in a local reference 
system attached to the barycentre are resulting from a set 
of consecutive transformations: 

(i) A rotation about the radial axis, ko, for aligning the 
first coordinate axis along the proper motion in the (io, jo) 
tangential plane (figure [TJ. The rotation angle j3 is derived 
from 

COS/3 = fl a */fl, 

sin/3 = fi s /fi. 

Let (ifo , vo , wo) be the new reference system, and let xo, 
%j)o and Co be the coordinates of the point in this system. 

(ii) The transformation of the (uo, vh, wb) reference sys- 
tem in (u,v,w), due to the displacement of the star. This 
transformation consists in a rotation about vo of an angle 
a, and in a correction of the angular coordinates when the 
distance of the barycentre, do, is changed in d (figurc[2]). a is 
the motion of the barycentre due to the proper motion, and, 
in a first order approximation, do/d — 1 — fi r t. Therefore, 
the coordinates in the (u, v, w) system are: 



= (1 — fi r t) 



where (xo,V>o,Co) refers to the (uo,vb,wo) frame. 

(iii) A rotation about the radial axis w must still be ap- 
plied in order to get the coordinates of the point in the most 
suitable reference frame for the forthcoming calculations. In 
Section 15.1.31 the motion of the binary is projected on a 
unique equatorial tangential plane, corresponding to a point 
arbitrarily chosen. This plane is called the 'average' equato- 
rial tangential plane hereafter; for simplicity, we chose, for 
reference frame, the local equatorial tangential frame of the 
beginning, (io,jo). Therefore, we need, for each barycentric 
position, an orientation of the (i, j) axes such that the pro- 
jections of these axes on the average equatorial tangential 
plane are parallel to (io, jo), respectively. Since the projec- 
tion of the motion of the barycentre is a straight line on the 
average plane, this condition is achieved by applying to the 
(u,v,w) frame a rotation with the angle —j3. The (i,j,k) 
frame thus obtained is called the 'local standard reference 
system' hereafter. 

In a first order calculation, sin a = fi t, cos a = 1, and 
the products between the angular motions, such as fi r t x 



fl a *t, are negligible. Therefore, the (x,y,z) coordinates in 
the (i,j, k) reference frame are given by the equations: 
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X = (1 — fl r t) X ~ {fl a *t]rad Z 
V = (1 — flrt) yo — [fJ-st]rad Zo 

Z = (1 - firt) Zo + [fl a *t] rad Xo + [flst}rad VO 



(1) 

(2) 
(3) 



where '[ ] ra d' means that the angle between the brackets is 
converted in radians (fi r t is also in radians, according to the 
definition of fi r above). 

An anonymous referee pointed out that the local standard 
reference system (i,j,k) is not a local equatorial reference 
system, contrarily to the (io,jo,ko) frame of the beginning 
of the calculation. The coordinates in the local equatorial 
reference system are not used in the derivation of the or- 
bital elements with the method of the average equatorial 
tangential plane, in Section [SJ but they may be relevant for 
high-velocity binary stars, when this method cannot be ap- 
plied. For that reason, they are derived hereafter. 

For getting the coordinates in the local equatorial refer- 
ence system of the barycentre, it would still be necessary to 
add a rotation about k with an angle equal to the rotation of 
the north, but in the opposite direction. If AO is the correc- 
tion to be added to the position angles measured inside the 
celestial sphere, as in the usual definition (figure[T}, then, the 
angle of the rotation to apply to the (i,j, k) reference frame 
for getting an equatorial frame is also AO. It comes from 
the b asic formulae of spherical trigonometry l|van de Kamrj 
Il967t ) that AO is given by the exact equation: 



tan A9 = 



sin a tan S tan 6 + (1 — cos a) sin 6 
1 — sin a tan 8 cos 6 — (1 — cos a) cos 2 6 



When a is small, the equation above may be simplified 
assuming (1 — cos a) = 0. If the binary is not close to a pole, 
it may also be assumed that <rtan<5cos# <C 1. Using the 
same approximations as above, the referee obtained: 

A6 = [fi a *i\radtan5 

Therefore, applying a rotation with a small A8 angle leads 
to the coordinates 

x\=x + [fi a ,t\ rad taxi5 yo 

yi = y — [fla*t]radt>aO-5 Xo 

in the local equatorial tangential plane of the barycentre. 



2.2 Application to an orbit 

2.2.1 The ascending node problem 

Two different sets of parameters are commonly used to en- 
tirely describe an orbit in space. Both are including the pe- 
riod, P, the eccentricity, e, and the periastron epoch, To. 
These parameters may be completed with the angular semi- 
major axis, a, and the Campbell's angles, which are: the 
inclination, i, the position angle of the ascending node, fi, 
and the so-called longitude of periastron, u>, which is mea- 
sured from the ascending node; the definitions of these an- 
gles may be found i n any double star handbook (Birmcndiik 
1960; lHeintzlll978h . Hereafter, the semi-major axis refers 
to the astrometric orbit of the primary component, which 
should not be confused with the relative 'visual' orbit of a 
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resolved binary system. The nodes are the intersections be- 
tween the orbit and the plane; the ascending node is 
the one where the star is moving away from the Sun. It is 
worth noticing that the positions of the nodes are derived 
from astrometric data, but that the ascending node is usu- 
ally not known when no RV measurement is available. For 
an astrometric orbit, Q and w refer then to the node with £1 
between and n. However, when the right ascending node 
is eventually found, it may be the other one, and Q and u 
are then changed in Q. + tt and u + n, respectively. 

The second set of parameters which are frequently used 
are the Thiele-Innes elements (A, B , F, G). These elements 
are functions of a and of the Campbell's angles, so it is 
possible to go fro m one system to the other and vice-versa 
l|Binnendiiklll960T ). However, changing Q and u in S~l + iv and 
u> + n doesn't change (A, B , F, G). Reciprocally, changing 
(A, B , F, G) in (a, i, Q,u>) leads again to the ambiguity in 
the choice of the node. 

The Thiele-Innes elements are completed with C and 
H , which are defined by the following equations: 



2.2.3 Fixing the ascending node from local perspective 



C — a sin i sin u> 
H — a sin i cos ui 



(4) 
(5) 



These two additional elements are not entirely redun- 
dant with the four former, since they indicate the true as- 
cending node, when it may be fixed. Otherwise, it is impos- 
sible to discriminate between (C,H) and (—C,—H). 



2.2.2 Correcting the Thiele-Innes elements for LP 

The six Thiele-Innes elements, (A, B ,C, F,G, H), are com- 
pletely fixing the orbit of a binary component in space. In the 
local reference system defined in Section 12.11 above, the di- 
rection cosines of the semi-major axis and of the semi-minor 
axis of the orbit are (B /a, A/a,C/a) and (G /a, F/a, H /a) 
respectively. The transformation of (A, B, F, G) by LP is 
obtained by applying equations ^ and <(2j) to these coordi- 
nates; one obtains: 



A = (1 — fX r t)Ao — [/J,st]rad Co 
B = (1 — fI r t)B() — [fla*t] ra d Co 
F = (1 - Hrt)F - \pgt\rad H 
G = (1 — fl r t)Go — [tJLa*t]rad Hq 



(6) 



It is worth noticing that the elements thus obtained are 
not the true Thiele-Innes elements, since the position angle 
of the line of nodes is not measured exactly from the local 
north, but from the j vector, ie from the north at t = 0. 

The position of a binary component is derived from 
the (A, B, F, G) elements with the classical formulae; in the 
(i,j,k) coordinate system defined in Section [2. II they are: 



x — B X + GY 
y = AX + FY 

with 

X(t) = cos E(t) - e 
Y(t) = VI -e 2 sin E(t) 

where E(t) is the eccentric anomaly at the epoch t. 



(7) 



(8) 



As explained in Section 12.2.11 when the actual ascending 
node is not known, the (Ao, Bo, Fo, Go) elements correspond 
to two sets (Co, Ho) : the first one with w£ [0, 7r[ and the sec- 
ond one with ui £ [tt, 2n[. However, thanks to the perspective 
effect, it is possible to discriminate between the two possibil- 
ities, using the following method: The (A, B, F, G) elements 
at each observation epoch are derived from equation ©, 
and the model positions of the binary component are de- 
rived twice, assuming the two u. The right ascending node 
will provide the solution with the smallest \ 2 ■ 

2.2.4 Size of the LP effect 

Our aim is to estimate of the size of the LP effect under the 
best conditions. For that purpose, we consider a binary at 
a distance of 5 pc. In order to get estimations applying for 
resolved binaries as well as for unresolved binaries, the mass 
ratio (mass of the secondary component divided by that of 
the primary one) q — 0.5 is assumed. For resolved binaries, 
the semi-major axis of the secondary orbit is then exactly 
twice that of the primary, and for unresolved binaries, the 
magnitude difference is large enough for neglecting the dif- 
ference between the position of the primary and that of the 
photocentre. The semi-major axis of the primary orbit is set 
to 1 AU; for a 1-solar mass primary and q = 0.5, this corre- 
sponds to a period of y/TE = 4.24 years. The assumed time 
span of the observations is 5 years, like the expected dura- 
tion of the Gaia mission. The orbit is then entirely covered 
by the observations. In order to choose a reasonable proper 
motion, a velocity of 25 km/s is assumed for each axis: this 
value is close to the standard de viation of the velocit y of 
late-type dwarfs, in any direction (Schaife rs et al.l [l982h . 

We consider now an edge-on orbit with the semi-major 
axis along the line of sight and the periastron between the 
barycentre and the sun. Therefore, Ao — Bo = 0, Co = —a; 
Assuming that the star is exactly at the opposite of the 
periastron when t =5 years, one obtains E = n, and equa- 
tion ([8]) leads to X — — (1 + e), Y = 0. Applying equa- 
tions and (0), the coordinates of the component are: 
x = —[fi a *t] rad a (1 + e), and y = ~[^st] rad a (1 + e), in- 
stead of when the perspect ive effect is absent. The offset 
in position in then x 2 + y 2 — 7.2 p.as for a nearly circular 
orbit. 

This estimation applies to a 1 AU semi-major axis at a 
distance d — 5 pc, but it may be adapted to other values. 
Since the Thiele-Innes elements are varying as aAu/d and 
the /x coordinates are varying as v km / a /d the LP effect is 
varying as Wfc m / S aAu/d 2 . Therefore, it is vanishing rapidly 
when distant stars are considered. 



3 THE LTT EFFECT AND THE ASCENDING 
NODE 

3.1 The LTT effect 

The light coming from a component orbiting around the 
barycentre of a binary system is arriving with a delay, 
tB + At, where tB is the delay of light coming from the 
barycentre, and where At is coming from the position of 
the component on the orbit. For that reason, At is called 
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the 'orbital delay' hereafter. It is negative when the com- 
ponent is closer to the Sun than the barycentre, and neg- 
ati ve otherwise. The analytic expr ession of At was derived 
bv lAnglada-Escude fc Torral l|2006n . For our purpose, a sim- 
ple calculation neglecting the relativistic corrections is suf- 
ficient, as it will be shown in Section 2] 

The radial angular coordinate with respect to the 
bary centre is derived from C and H with the following equa- 
tion (|Couteaulll978r ): 



C(cos£-e) + Hy/l - e 2 sinE 
Therefore, the orbital delay is 



At 



(9) 



(10) 



where c is expressed in AU per time unit, and where zu is 
in the same unit as z, i.e. as C and H. An astrometric mea- 
surement obtained at the epoch t then provides the position 
of the component at the epoch t — At. This induces a shift 
of the position of the star which is related to the orbital mo- 
tion, but also to the proper motion of the barycentre (see 
the Fig. 2 in Anglada-Escude & Torra). 

As for LP, the choice of the true ascending node is fea- 
sible trying the two possibilities. Shifting from one node to 
the other consists in changing the sign of z, and then of At. 
Again, the ascending node should provide the best fit when 
this effect is visible. 



3.2 Size of the LTT effect 

For comparing the offset of position due to the orbital de- 
lay to the LP effect evaluated in Section 12.2.41 the same 
hypotheses are assumed. With an edge-on orbit, the orbital 
delay is then as large as 8.3 minutes at maximum. The cor- 
rections of position coming from the orbital motion and from 
the proper motion are evaluated hereafter. 



3.2.1 The orbital motion 

The orbital velocity is given by the equation: 



VAU/y 



2tt 



M q 



3/2 



oau 1 



where M is the mass of the primary star, in solar unit; 
A4 = 1 is assumed. With q = 0.5, and a period of 4.24 y, 
v Au/y — 1-5- For 8-3 min, the correction is then 4.7 u.as for 
the primary component at a 5 pc distance. When another 
distance or another semi-major axis (i.e. another period) are 
considered, since At is proportional to oau, this correction 
is varying as y/a~Au / d pc . Again, it is also twice larger for the 
secondary component than for the primary. 



3.2.2 The proper motion 

Since the primary component is observed with an orbital 
delay At, the position of the barycentre which is used as 
origin of the orbital motion is no more corresponding to the 
epoch t, but to t — At. The apparent orbit is then distorted, 
since At is varying in relation with the orbital motion. 

For consistency with Section ^. 2. 41 a tangential velocity 
of 25\/2 = 35 km/s is assumed. Therefore, the correction of 



the position for a 8.3 min delay is 24 |i.as for a distance of 
5 pc. This term is varying as v km / s x a,Au/dp C - 

It is worth noticing that the parallactic motion is not 
affected by At, since it is due to the position of the Earth 
with respect to the Sun at the time of the observation. 



4 LOCAL RELATIVISTIC EFFECTS 

In Sections [2] and [3] above, the relativistic effects due to the 
displacement of the binary with respect to the barycentre of 
the solar system were ignored. In this Section, it is verified 
that this approximation doesn't matter to the problem we 
are dealing with. 

The linear motion of the barycentre of the binary has 
two consequences: the lengths along the direction of the dis- 
placement are shortened, and time is dilated. The coordi- 
nate of the primary co mponent al ong the velocity vector, v, 
is then multiplied by \J\ — v 2 /c 2 . Therefore, with the same 
hypotheses as above (spatial velocity of 25\/3 = 43 km/s, 
orbital radius of 1 AU, distance of 5 pc) we get a change in 
position that may be as large as 2 10 -3 \i&s. This is much 
smaller than the LP or LTT effects. 

The effect of time dilatation is as if, for the observer, 
the period of the binary is increased by — v 2 /c 2 . This 

linear modification of time has no effect on our calculations. 
Anyway, it induces an amount of only 1.4 s for a 4.2 year 
period. 

We may also consider the deflection of light by the mass 
of the secondary component. For a distant binar y object, we 
use t he approximate formula of light bending (jStraumannl 
1 19841): 



4G 



M 

' c 2 b 



where G is the gravitational constant, M is the mass of the 
deflecting star, and b is the 'impact parameter' of the ray, 
measured from the centre of the deflecting star. Therefore, 
in the reference system defined in Section [5] the deflection 
of the light emitted by the primary component by the sec- 
ondary component is, approximately: 



= 4G 



M.2 *1|2 



where zi\2 is the radial angular coordinate of the observed 
component (labelled '1') in a local reference system with 
component 2 as origin; therefore, deflection occurs only when 
Zi|2 > 0. With the same hypotheses as before, we have M2 = 
0.5 solar mass, z^ 2 = 3 AU / 5 pc = 0.6 arcsec. The light 
deflection may then be as large as around 1.3 (i.as when 
b — 2 solar radii. This is not negligible, but corresponds to a 
rather rare situation, i.e. an orbit almost perfectly edge-on. 
Since 6<f> is varying as 1/b, it vanishes very fast when b is 
increasing from a very small value. When the median value 
of 7r/3 radians is assumed for the orbital inclination, Scj> w 
0.007 u.as only. In conclusion, light bending may be neglected 
in a first approximation, although it may be important for 
some systems. 
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5 LP AND LTT EFFECTS IN SIMULATED 
UNRESOLVED BINARIES 

It comes from Section [2T2T41 arid I3T21 that LP and LTT effects 
could significantly alter the position of the photocentre of 
an unresolved binary when the accuracy is around 10 u.as. 
However, simulations are necessary to evaluate when the 
signatures of these effects should be visible in reality. 

In this section, we attend to answer the two questions 
addressed in the introduction: 1) how does the distribution 
of the GOF look when LP and LTT are neglected?, and 2) 
When is it possible to select the actual ascending node in 
practice? 

5.1 The simulated binaries and their observations 

5.1.1 Properties of the virtual systems 

Several hypotheses in 12.2.41 are assumed again for gener- 
ating virtual binaries: the mass of the primary component 
is 1 solar mass and that of the invisible companion is 0.5 
solar mass. The period is 4.24 years, so that the semi-major 
axis of the astrometric orbit is exactly 1 AU; the distance is 
5 pc. In addition, each coordinate of the barycentric velocity 
obeys a normal distribution with the standard deviation of 
25 km/s, the positions of the stars are randomly distributed 
on the sky, and the orbits are not circular, but have eccen - 
tricities obeying the distribution of lHalbwachs et ail (2005). 
The orientations of the orbits are purely at random: the in- 
clinations obey a sinus distribution, and the Q and ui angles 
obey constant distributions. 

5.1.2 Properties of the observations 

In order to reproduce the data coming from a rotating satel- 
lite or an interferometer, we consider one-dimensional astro- 
metric coordinates ('abscissae' hereafter) measured along an 
axis with a varying orientation. For a given binary system, 
the origin of the abscissa corresponds to a fixed position close 
to the position of the barycentre at t — (Fig. |4]). Virtual 
observations are produced assuming a rough approximation 
of the Gaia scanning law, consisting in around 16 visibility 
periods (VP, hereafter) in 5 years, separated by a duration a 
bit longer than 100 days. A constant orientation is randomly 
fixed for each VP, which contains 2 or 3 pairs of scans. The 
total number of measurements is 83 on average. Each pair 
corresponds to a complete rotation of the satellite in 6 hours. 
The second scan of the pair is coming 106 min after the first 
one. For any scan, an abscissa is derived, and an error is gen- 
erated from a normal distribution with a constant standard 
deviation, o~ w . 

5.1.3 Deriving the abscissae from the parameters of the 
stars 

The easiest way to derive the parameter of a star is to con- 
sider its motion in the average equatorial tangential plane 
defined in Section 12.11 When the distance as between the 
origin and the barycentre of the star is small, the curvature 
of the celestial sphere may be neglected, and it is assumed 
that tan erg = as - The coordinates of a binary component in 
the local standard tangential plane, seen in Section [2~T1 




to sun 



Figure 3. Cross section of the celestial sphere, from the local 
standard tangential plane to the average equatorial tangen- 

tial plane (i a ,j s ). B is the barycentre of the binary as seen on the 
celestial sphere, and P is a component or a peculiar point of the 
orbit. The angular distance between B and P is the arc Bp. It is 
nearly equal to ||Bp'|| rj ||bp" ||, where p' is the projection of P on 
the plane, and ||bp"|| is the projection of Bp on the average 

equatorial tangential plane. 

may then be directly transposed to the average equatorial 
tangential plane, neglecting the tilt with the angle as be- 
tween the two planes (see Figure [3}; they are just added to 
the coordinates of the barycentre. 

However, we must verify that the curvature of the celes- 
tial sphere may really be neglected. Assuming a 5 pc distance 
and a velocity as large as 2a on each axis in the tangential 
plane, one obtain, /i — 3.0 arcsec yr _1 . For half the duration 
of a 5 yr mission, this implies a displacement of 7.5 arcsec. 
The difference between the angular distance and its projec- 
tion on the tangent plane is then only 0.003 pas. It is then 
quite possible to use the plane approximation, even when 
astrometric errors as small as 0.1 pas are considered. 

The calculation of the abscissa of the star for the ob- 
servation epoch t is performed in several steps, which are 
roughly summarised in Fig. U Hereafter, LP and LTT af- 
fects are introduced in the calculation; about the orbital 
delay, since it is much smaller than the period of the binary, 
a first order approximation is used: The eccentric anomaly 
E{t) is computed, and the z-coordinate and the orbital delay 
At are derived, from equations Q and (|10|l . The eccentric 
anomaly corrected for LTT, E{t— At), is then computed, and 
the x and y coordinates in an equatorial rectangular frame 
are derived from equation (0. The contribution of the or- 
bital motion to the abscissa of the star, w , is the projection 
of the (x,y) coordinates on the abscissa axis. When <f> is the 
position angle of the axis, w is given by the expression: 

w — x sin <j> + y cos <j) (11) 

The abscissa of the star is then derived by still adding 
to w the position of the barycentre, wb ■ The perspective ef- 
fects concerning single stars are then taken into account: the 
parallax is varying due to the radial motion, and the proper 
motio n is affected by perspective acceleration (|Dravins et al.l 
Il999l ). Therefore, the abscissa of the barycentre, wb, is given 
by the equation: 
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Figure 4. The abscissa of a binary component in the average 
equatorial tangential reference plane. It is the projection of the 
heliocentric position of the barycentre at t = (indicated by the 
cross), plus contributions of the proper proper motion, of the par- 
allax, and of the orbital motion of the star around B, the position 
of the barycentre at the t epoch. The origin of the rectangular 
equatorial coordinates is an a priori rough estimation of the av- 
erage position of the star. 



wb ~ Aa*o sin <f> + ASo cos 4> + (1 — fi r x (t — At)) x 

(ro-o/ro + Ma*o x (* - At) sin (j) + (12) 
fjbs x (t — At) cos 4>) 

where Aa*o, A<5o, vjo, /i Q * , and fis are the position, paral- 
lax and proper motion of the barycentre at the epoch t = 0, 
respectively. / ra is the partial derivative of w with respect 
to the actual parallax at the observation epoch t. It depends 
on the coordinates of the binary and on the position of the 
Earth on its orbit. The (1 — fj, r X (t — At)) coefficient is due 
to barycentric perspective. The introduction of barycentric 
perspective in the simulated wb is necessary for being real- 
istic, and it is very relevant in practice: in the next section, 
it will make possible the derivation of the RV with the other 
parameters; therefore, it contributes efficiently to the cor- 
rection of the LP effect when the RV is not measured spec- 
troscopically. It is worth noticing that the effect of the epoch 
correction in this term is negligible, and that it is possible 
to simply use (1 — /x r t). 

Developing w with equations (JSJ) to (J7]), one obtains 
the sum w + wb '■ 

w(t) = Aa*o sin cj> + A<5o cos (j> 

+ roo/ra (1 — Urt) 

+ [i a * (1 - [t, T t - [z(t - At)] rad ) (t - At) sin <p 

+ ^S a (1 — fMrt - [z(t— At)]rad) (t — At) COS <f) 

+ A (1 - Urt) X(t - At) cos cj> (13) 
+ B (1 - Urt) X(t - At) sin <j> 
+ F Q (1 - fi r t) Y(t - At) cos4> 
+ Go (1 - Urt) Y(t - At) sin <t> 



which is the abscissa of the photocentre when LP and LTT 
are both taken into account. It is worth noticing that equa- 
tions (|lip to (|13|) are not specific to an observation mode, 
and apply to any coordinate axis. When 2-D measurements 
are considered, the ordinates obey the same formulae as the 



abscissae; only one parameter is changing: the position angle 
of the coordinate axis, <j>. 

5.2 Calculation of the solutions with local 
perspective and LTT 

Thirteen parameters are necessary to entirely fit the mo- 
tion of the primary component of a binary system, plus the 
indication of the ascending node. In the usual order, these 
parameters are: 

• The 6 parameters of the barycentre of the binary at the 
reference epoch t = 0. These parameters are the five ones 
listed in Section [5.1.31 above (i.e. Aq.o, ASo, ^o, Ma* > fJ-s ), 
and the RV, v r . Contrarily to the parallax and the proper 
motion, the RV is assumed to be constant. This parameter 
is optional in the calculation of the solution, since it may be 
obtained independently from the astrometric data; for in- 
stance, it will be derived from spectroscopic observations in 
the Gaia project. Therefore, two different cases are consid- 
ered hereafter: the case with known RV, and the case with- 
out velocity measurement, for which v r is included among 
the unknowns of the astrometric solution. 

• The Thiele-Innes elements at the reference epoch t = 
0, which are (Ao, Bo, Fo, Go). Two sets of the elements Co 
and Ho may correspond to these 4 elements, as indicated 
in equations and For each of them, the elements 
(A, B, F, G) are given by equation ([6} for any epoch (t — At). 

• The period, P, the eccentricity, e, and the epoch of 
periastron, To- These elements are necessary to derive the 
eccentric anomalies, E(t — At), and therefore the 3-D coor- 
dinates of the binary components, (x,y,z). 

These parameters are derived by x 2 minimisation, us- 
ing the Lev enberg-Marquardt algorithm and the routines 
provided bv lPress et al.l (|1992l ). 

Each time a solution is calcula ted, its gaussianized GOF 
is derived, using the F% estimator (|Stuart fc OrJl994lESAl 
[l997|): 



n I 9n 



(14) 



where n is the number of degrees of freedom. 

For each virtual binary and its set of 'measured' abscis- 
sae, the astrometric and orbital elements are derived three 
times, always taking into account the barycentric perspec- 
tive effects: a first solution is derived from a simple orbital 
model, ignoring the local effects related to the z-dimension. 
LP and LTT are taken into account in the two other solu- 
tions, which are obtained assuming the right ascending node, 
and the other one, respectively. The quantity AFb, defined 
as the difference: 'GOF of the false ascending node solution' 
minus 'GOF of the right ascending node solution' is then 
derived. 



5.3 The rejection risk when local perspective and 
LTT are neglected 

Our aim is to estimate a limit of astrometric errors beyond 
which neglecting LP and LTT would increase the rate of 
rejected solutions. For that purpose, the solutions derived 
without LP and LTT are considered. The distributions of 




Figure 5. The GOF of the solutions derived for the simulated 
binaries, when LP and LTT are neglected. The actual radial ve- 
locities are assumed for the barycentres of the systems. The dis- 
tribution clearly deviates from the normal Af(0, 1) distribution 
when the error of the astrometric measurements is less than 1 uas. 



Figure 7. Distribution of the differences between the F2 esti- 
mators of the GOF when the actual radial velocity is assumed. 
AF2 < means that the solution with the ascending node in the 
right place leads to a x 2 > and therefore a GOF, larger than the 
solution with the node at the opposite place. 




Figure 6. Same as Fig. [5] but including the radial velocity among 
the unknowns of the solutions. The distribution is no more a 
normal Af(0, 1) distribution when the error of the astrometric 
measurements is less than 0.5 uas. 



the GOF when the actual RV is assumed is plotted in Fig. [5] 
Fig. HJ] refers to the solutions obtained including the RV 
among the unknowns. 

On both figures, it appears that a 10 uas error doesn't 
alter the distribution of the F2 estimator, which just looks 
like the normal distribution N(0, 1). A significant excess of 
systems with F2 > 3 appears only when cr w is less than 1 uas 
when the RV is known, and when <r w < 0.5 u.as otherwise. 

The errors leading to visible LP or LTT effects are 
rather small when compared to the estimations in Sec- 
tion [2T2T4] and [372] The reason is that the offsets of positions 
due to these effects are partly compensated by altering the 
orbital elements. In fact, the modification of the orbital el- 
ements is so efficient that the \ 2 °f the solution doesn't 
change significantly when it is accounted for the LTT effect. 
This appears when simulations are performed assuming null 
velocities for the barycentres of the systems, in order to sup- 
press the LP effect. The GOF of the solutions are then al- 
ways obeying the same normal distribution, no matter such 
effect is taken into account, and no matter the precise choice 




Figure 8. Same as Fig[7] but when then the RV is derived from 
the astrometric data. 



of the ascending node. On the other hand, the errors of the 
derived eccentricities and of the Campbell's angles are much 
larger when the wrong node is assumed in place of the right 
one: for instance, for the eccentricities, when <r w = 1 uas, 
the ratio is already more than 2. Therefore, taking LTT into 
account improves the quality of the solution, but it doesn't 
help in the selection of the ascending node. It is worth notic- 
ing that this conclusion applies to unresolved binaries, and 
that the LTT effect is more visible in resolved binaries, as 
explained in Section [6] 



The conditions for selecting the ascending 
node 



5.4 



The differences AF2 defined in Section [5T2l were derived from 
the solutions taking LP and LTT into account, and their 
distributions are plotted in Fig.[7]and|8] It comes from these 
plots that, for some systems, F2 may be used to select the 
ascending node even when the astrometric error is too large 
for providing an excess of large GOF in Fig. [5] and [6] For 
instance, when a w = 1 uas, the distribution of AF2 is clearly 
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Figure 9. Standard deviation of the radial velocity error, as a 
function of the error of the astrometric measurement. 



asymmetric, and |Ai<2| > 0.5 means that the true ascending 
node is very probably the one with the smallest F2. 

5.5 The radial velocity derived from perspective 

Not surprisingly, the solutions derived assuming the actual 
RV of the barycentre are more efficient in the node selection 
than that including the calculation of v r . Since the actual 
RV can never be measured in practice, the results above 
can only refer to RV 'accurate enough'; but what accuracy 
may still be considered as 'enough' ? In order to answer this 
question, we consider the RV errors when v r is derived from 
the astrometric measurements. The standard deviation of 
these errors, a^RV, is plotted in Fig. H as a function of 
the astrometric error, a w . An average RV measured with an 
uncertainty much smaller than uarv may be considered as 
equivalent to the actual one. 

In the Gaia project, it is expected to get the average RV 
of the brightest systems, with an accuracy of a few km/s. 
This is much less than the RV error when v r is derived from 
the astrometric measurements, since the astrometric error of 
each Gaia measurement will be much larger than the upper 
limit of the range of a w considered in Fig. [§] Therefore, 
the possibility to find the true ascending node with Gaia 
could be obtained through the solutions assuming the actual 
RV. Unfortunately, since the error of the astrometric Gaia 
measurements at epoch (not to be confused with the average 
astrometric error) will be much larger than 10 u.as, and it 
is obvious that LP and LTT may both be ignored in the 
preparation of the Gaia catalogue. 



6 APPLICATION TO RESOLVED 
ASTROMETRIC BINARIES 

6.1 Calculation of the astrometric solution 

The resolved astrometric binaries are less frequent than the 
unresolved ones, but they are also more relevant, since the 
masses of the components are then directly obtained from 
the astrometric solution. The derivation of the parameters 
is not very different from the calculation seen above: the ab- 
scissae of the primary component are again given by equa- 
tion (|13[) . and the method described in Section [5.21 may be 



applied with the parameters of the primary component in 
place of that of the photocentre. 

The abscissae of the secondary component are slightly 
different. It is necessary to add a parameter to that of the 
primary, in order to derive the Thiele-Innes elements of the 
secondary. The simplest method consists in adding the sec- 
ondary semi-major axis, 02- The six Thiele-Innes elements 
of the secondary may then be derived from that of the pri- 
mary, since A2 = — 02-Ai/ai, B2 = — a2-Bi/ai, and so on, 
until H2 = — a2-Hi/ «i (a method for deriving ai from Ai, 
Bi, Fi and Gi is in lBinnendiik|[l96rl . Therefore, X, Y, z 
and At may be derived for the secondary using the same 
formulae as for the primary. 

When the parameters of the binary have been com- 
puted, the masses of the components are derived from 




(12 



Oi 



(01+02) 
(01+02) 



(15) 



where the masses are in solar units when 01 and a2 axe in the 
same unit as zu. Po is the period at rest, expressed in years; 
it is derived from the apparent period, after correction for 
the Doppler effect: P = P/(l + v r /c). 



6.2 Application to 2MASS J07464256 +2000321 

6.2.1 The LP+LTT solution compared to the 'no local 
effect' solution 

Instead of repeating the calculations performed for unre- 
solved binaries, it is more appealing to consider a real nearby 
binary, and to search when LP and LTT must be taken into 
account in order to get masses as accurate as possible. The 
sample of nearby systems for which the orbital elements are 
already known is rather scarce. 2MASS J07464256 +2000321 
roughly fulfils the requirements: it is a reasonably nearby 
system (12 pc distant), with a period long enough for hav- 
ing a wide angular separation between the component (semi- 
major axis of 0.2 arcsec), but short enough (10.5 y) for per- 
mitting observations covering a complete orbit. Neverthe- 
less, the observation frequency used in Section [5] is adapted 
hereafter in order to keep the same number of observations 
although the time-span is extended up to 12 years. The pa- 
rameters of 2MASS J07464256 +20003 21 are summa r ised i n 
Table 1. The radial velocity comes from lBailer-JonesI d2004h 
the m asses and the orbital parameters are from lBouv et al.l 
(2004); the other data are from Simbad. 

Synthetic observations of 2MASS J07464256 +2000321 
were generated, assuming a constant astrometric error for 
each run. The difference between the beginning of the ob- 
servations and the periastron epoch was randomly varying 
between and P. 

Three solutions were computed for each set of virtual 
observations: 

(i) the solution with the ascending node in the right place, 
which is called the 'LP+LTT solution' hereafter; 

(ii) the solution with the ascending node at the opposite 
place; this solution was used to derive the difference AF2; 

(iii) a solution ignoring LP and LTT effects, which is re- 
ferred as the 'no local effect solution'. 

For the solutions (i) and (ii), the errors of the masses, 
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Table 2. The node solution and the 'no local effect' solution for 2MASS J07464256 +2000321, for various 
astromctric errors. 



a w LP+LTT solution no local effect solution 



(Mas) 


o(AMi) 


v(AM 2 ) 


AF 2 > 


(AMi) 


a(AMi) 


(AM 2 ) 


a{AM 2 ) 


(F2) 




(Me) 


(M e ) 


(per cent) 


(M @ ) 


(Mq) 


(Mq) 


(M @ ) 




1 


7.7 10~ 7 


6.1 10~ 7 


99.5 


-11 10~ 7 


22 10~ 7 


-7.2 10~ 7 


17 10" 7 


6.8 


2 


1.5 10" 6 


1.2 10~ 6 


92 


-1.1 10~ 6 


2.6 10~ 6 


-0.7 10~ 6 


2.0 10" 6 


2.1 


5 


3.8 10" 6 


3.1 10" 6 


72 


-1.1 10~ 6 


4.4 10" 6 


-0.7 lO" 6 


3.5 10" 6 


0.37 


20 


7.7 10~ 6 


6.1 lO" 6 


61 


-1.1 10" 6 


8.0 10~ 6 


-0.7 10- 6 


6.4 10" 6 


0.10 


50 


3.8 10~ 5 


3.1 10~ 5 


53 


-0.1 10" 5 


3.8 10- s 


-0.1 10~ s 


3.1 10" 5 


0.01 



Table 1. The astrometric elements of 2MASS J07464256 
+2000321, as assumed in the simulations. 



2MASS J07464256 +2000321 





81.9 mas 


Ha* 


-368 mas/y 




-39 mas/y 


V r 


54.1 km/s 


I' 


10.5432 y 


e 


0.41 


i 


141.6 dcg 


Q 


20.7 deg^. 


iO 


350.6 deg 


Mi 


0.085 Mq 


M 2 


0.066 Mq 



a f2 is known only modulo n. 



A All and AA4 2 were computed, as well as their standard 
deviations, o. The errors are positive when the solution pro- 
vides masses larger than the actual ones, and negative oth- 
erwise. The F 2 estimator of the GOF was also computed for 
each of the three solutions, and the difference of F 2 between 
solution (ii) and solution (i) was derive; as in the previous 
section, this difference is positive when the x 2 of the solu- 
tion with the right ascending node is smaller than that of 
the other node. The most relevant statistics are gathered in 
Table 2. 

It comes from these results that F 2 is much more sensi- 
tive to the local effects than for an unresolved binary with 
the same semi-major axis. The reason lies in the compensa- 
tion of the LTT effect by alteration of the orbital parameters, 
as explained in Section 15.31 For resolved binaries, the sign 
of At is not the same for both components; therefore, ig- 
noring LTT results in shifting the position of the barycentre 
out from the line joining the components; no modification 
of the orbital parameters can reproduce such alteration. 

Another important result visible in Table 2 is the sys- 
tematic error when the local effects are neglected. The stan- 
dard deviation of the mass errors are becoming similar when 
the astrometric error, a w , is around 20 u.as, but the average 
error is still 14 per cent of the error of the LP+LTT solution. 
The LP and LTT effect are really negligible only when a w 
is around 50 |i.as. 



Table 3. The mass errors and the GOF coming from the degen- 
erate LP solution for 2MASS J07464256 +2000321. 





(AMi) 


<t(AA4i) 


(AM 2 ) 


<t(AM 2 ) 


<F2> 


(uas) 


(Mq) 


(Me) 


(Mq) 


(Mq) 




1 


-2.4 10~ 7 


9.4 10" 7 


-0.7 10~ 7 


7.4 10" 7 


0.48 


2 


-2.4 10~ 7 


1.6 10" 6 


-0.7 10~ 7 


1.3 10~ 6 


0.13 


5 


-2.4 10" 7 


3.9 10" 6 


-0.7 10~ 7 


3.1 10~ 6 


0.03 



6.2.2 The 'degenerate LP solution' 

Accurate masses may be computed with a LP+LTT solu- 
tion, but only when the position angle of the ascending 
node is known. At present, the ascending node of 2MASS 
J07464256 +2000321 is still not fixed despite of several ra- 
dial velocity measurements, because t he spectrum of th is 
star is a blend of the two components ( Bai ler- Jonesll2004h . 
The comparison of the GOF of the LP+LTT solutions de- 
rived from the two nodes may lead to the selection of the 
ascending one, but it comes from Table 2 that this method 
would be reliable only when a w is less than 2 ^ias. There- 
fore, another calculation method must be used between 2 
and 20 (i.as. This method takes into account a part of the 
LP effect, but not the effects related to the position of the 
true ascending node. For that reason, its solution is called 
the 'degenerate LP solution' hereafter. 

The degenerate LP solution takes into account the vari- 
ation of the size of the orbit (i.e. of Ai t2 , B\ t2 , Fi t2 and Gi, 2 ) 
when the distance is changing. The calculation consists in 
ignoring z and At in equation (|13|1 . but not the (1 — fi r t) co- 
efficient of the Thiele-Innes elements. The results of the de- 
generate LP solution are presented in Table 3. As expected, 
they are quite acceptable for a w larger than 2 u.as, i.e. when 
the ascending node cannot be determined. 



7 SUMMARY AND CONCLUSION 

We have study the impact of two effects which have been 
neglected up to now in astrometric studies of binary sys- 
tems: the local perspective, and the orbital light -travel time. 
These effects are related to the orientation of the orbit in 
space, and taking them into account may lead to the posi- 
tion of the ascending node of the orbit, i.e. to the orientation 
of the spin vector. The statistical distribution of the spin ori- 
entations is a clue in the study of bina ry formation. It was 
derived and discussed in the past by iDommangetl (|l988l ) 
and, more recently by iGlebockJ <|2000l) . The global distri- 
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bution of the spin looks isotropic, but Glebocki concluded 
from a sample of 252 orbits that subgroups seem to have 
asymmetric distribution. An extension of his sample would 
allow to re-consider this question. 

We have found that neglecting LTT doesn't affect the 
X 2 of the solution for unresolved binaries, since it is com- 
pensated by a bias in the orbital elements, essentially the 
eccentricity and the orientation of the orbit. The reason is 
that LTT doesn't change the shape of the apparent orbit; 
only the observation epochs are affected, in relation with the 
orbital phase. Therefore, neglecting LTT can't lead to dis- 
carding the binary model. On the other hand, taking it into 
account for stars with distances d = 5 pc would significantly 
improve the accuracy of the orbital elements when the as- 
trometric accuracy a m is a few \ias. For other distances, the 
minimum a w is varying as the inverse of d. 

The ascending node of unresolved astrometric binary 
orbits may be found from the LP effect. However, for that 
purpose, astrometric measurements with accuracy around 
1 u.as and beyond are required. When the barycentric radial 
velocities are known, and therefore taken into account in the 
correction of perspective acceleration, LP must be taken into 
account when the errors of the astrometric measurements 
are less than the 1 [ias limit. The binary star model could be 
erroneously discarded otherwise, since the GOF of the solu- 
tion would look excessive. When the radial velocities are not 
available, they are derived from the barycentric perspective 
effects; however the corrections for LP are then necessary 
only when the errors of the astrometric measurements are 
less than 0.5 u.as. As for LTT, the results depending on LP 
apply to binaries at a distance d — 5 pc, with velocities typ- 
ical of population I late-type dwarfs, and with semi-major 
axes of 1 AU. For other parameters, they are unchanged 
when the astrometric errors are varying in the same propor- 
tions as the velocities and the semi-major axes, and in the 
inverse proportion of d 2 . 

The resolved binaries were examined through the exam- 
ple of 2MASS J07464256 +2000321, a brown dwarf system 
for which all the parameters are already known, the ascend- 
ing node excepted. The LTT effect is not compensated by 
alteration of the orbital elements as for unresolved binaries, 
since it plays in opposite directions for both components. It 
appeared that neglecting LP and LTT would lead to erro- 
neous masses when the errors of the astrometric measure- 
ments are around 20 [ias. However, when the astrometric 
error is larger than 1 or 2 u.as acceptable evaluations of the 
masses may be obtained from a 'degenerate LP solution'. 
This calculation doesn't require the knowledge of the as- 
cending node since it consists in correcting the semi-major 
axis for distance variations. However, when the astrometric 
error is around 1 (ias, the degenerate LP solution is no more 
acceptable, and only a calculation entirely taking into ac- 
count LP and LTT may lead to masses as accurate as per- 
mitted, through the determination of the ascending node. 
Since these results apply to a specific binary at a distance of 
12 pc and with very light components, it is highly probable 
that a few binaries in the forthcoming Gaia mission will re- 
quire a degenerate LP solution for getting accurate masses. 
On the other hand, complete LP+LTT solutions will be nec- 
essary only when measurements with errors around 1 nas 
will be available for nearby binaries, possibly thanks to the 
SIM project. 
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